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SECTION  I 
INTRODUCTION 

> 

The  use  of  a  large,  slow  energy  storage  device  to  store  energy  in  the  magnetic 
i  field  of  an  inductor  and  subsequently  transferring  a  significant  fraction  of  the  stored 

energy  to  a  parallel  load  by  means  of  an  opening  and  a  closing  switch  has  been  sug¬ 
gested  by  several  authors  (references  1-4)  and  is  illustrated  in  figure  1.  The  energy 
delivery  time  to  the  parallel  load  can  be  significantly  shorter  than  the  charging  time 
for  the  storage  inductor,  resulting  in  energy  pulse  compressions  of  factors  of  5  to  10. 

For  a  matched  load  (load  inductance  equal  to  storage  inductance),  a  maximum  of  25% 
of  the  initial  energy  in  the  storage  inductor  can  be  transferred  to  the  load  (reference  3) . 

If  the  load  is  dissipative  (for  example,  a  SHIVA  electromagnetic  implosion,  reference 
5),  a  substantially  higher  fraction  of  the  stored  energy  could  be  transferred.  A  fuse 
could  serve  for  the  opening  switch,  and  it  would  be  required  to  interrupt  currents  of 
many  megamps  and  to  hold  off  voltages  of  hundreds  of  kilovolts. 

Using  a  very  large  capacitor  bank,  homopolar  generator  (reference  6),  or 
explosive  magnetic  flux  compression  generator  (reference  7)  with  one  or  more  storage 
inductor/  opening  switch  pulse  sharpening  stages,  extremely  high  energy  (~  107  to 
10*  j)  and  high  power  (~  1014  to  10,3w)  pulses  may  be  possible.  Fuses  for  such 
pulsed  power  devices  would  interrupt  ~30  to  300  MA,  hold  off  ~300  kV  to  3  MV, 
and  absorb  at  least  10%  of  the  stored  energy  without  restriking. 

The  purpose  of  this  work  was  to  extend  the  fuse  performance  of  copper  foils 
achieved  by  Di  Marco  (reference  4,  800-kA,  80-kV,  25-kj)  to  a  2-MA,  250-kV,  200- 
kj  operation.  Also,  information  on  the  scaling  of  fuse  behavior  and  resistivity 
(references  4  and  8)  with  energy  was  sought. 

Some  useful  criteria  for  estimating  optimum  fuse  size  for  a  given  circuit  are 

those  given  by  Maisonnier  (reference  3).  The  first  criterion  is  that  the  energy 

2 

deposited  in  the  fuse  by  I  R  heating  should  be  just  enough  to  initiate  fuse  vaporization 
at  maximum  current  when  the  charge  left  stored  in  the  capacitor  bank  is  zero.  The 
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Figure  1 .  Schematic  Diagram  of  Inductive  Energy-Storage 
Circuit. 
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second  is  that  the  fuse  mass  should  be  such  that  the  energy  which  must  be  dissipated 
in  the  fuse  is  just  enough  to  vaporize  it.  These  criteria  result  in  conditions  on  the 
fuse  length  and  cross  section  in  terms  of  the  bank  capacitance,  total  inductance, 
initial  charge  on  the  bank,  and  the  fuse  material  properties.  These  criteria  were 
used  as  a  guide  for  a  parameter  survey  in  this  work. 


SECTION  n 


EXPERIMENTAL  APPARATUS 


The  experiments  consisted  of  discharging  a  capacitor  bank  into  a  fuse  package 
and  measuring  the  discharge  current  and  the  voltage  across  the  fuse.  The  capacitor 
bank  consists  of  4  modules  of  16  capacitors  each.  Two  modules  are  connected  in 
parallel  to  form  an  arm,  and  each  arm  discharges  through  four  pressurized  rail  gap 
switches  in  parallel.  The  outputs  of  both  arms  are  connected  in  parallel  to  the  fuse. 

The  measured  capacitance  of  the  bank  is  158  +  3  //F,  and  it  is  charged  to  50  kV 
to  give  an  initial  stored  energy  of  nearly  200  kj.  Short  circuit  shots  were  fired  in 
order  to  measure  the  inductance  of  the  bank  and  transmission  lines.  Excluding  the 
fuse,  the  system  inductance  is  46  +  2  nH. 

Figure  2  is  a  cross-sectional  view  of  a  typical  fuse  package.  A  rectangular 
piece  of  copper  foil  is  clamped  to  the  transmission  lines  at  opposite  ends.  The  foil 
is  folded  around  che  insulator,  but  it  is  not  creased.  The  two  halves  of  the  foil  are 
parallel  and  separated  by  1.5  cm.  A  polyethylene  envelope  is  taped  around  the  foil, 
but  the  envelope  does  not  touch  the  foil.  Glass  beads  (sand)  are  poured  into  the 
envelope  to  surround  the  foil.  Care  is  taken  to  ensure  that  the  beads  fill  the  space 
between  the  insulator  and  the  foil  on  both  sides. 

"Blast-O-Lite"  BT-12  glass  beads,  which  are  between  62  and  105//  in  diameter, 
were  used.  It  was  found  that  the  finer  BT-13  beads  (44  to  88  //)  showed  no  difference 
in  performance,  but  that  the  coarser  BT-8  beads  (149  to  210  jj)  resulted  in  signifi¬ 
cantly  longer  turn-off  times.  Since  the  BT-13  beads  are  much  more  difficult  to 
handle,  the  BT-12  beads  were  routinely  used. 

All  of  the  fuses  were  made  from  one  roll  of  rolled  copper  sheet.  An  electrolytic 
tough-pitch  copper  (alloy  110)  that  was  nominally  1-mil  thick  was  used.  Samples  of 
the  foil  were  weighed,  and  it  was  found  to  be  21.4  mg/cm2  thick.  Therefore,  the 
equivalent  thickness  of  the  foil  is  23.9  H. 
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Seventeen  different  fuse  sizes  were  fired  at  the  200-kj  level.  Sizes  were  deter¬ 
mined  by  varying  the  length  and  width.  The  ranges  were  up  to  +40%  from  a  reference 
size,  which  was  determined  from  Maisonnier's  equations  (reference  3).  These 
equations  employ  three  factors:  lq  ,  k2,  and  a,  which  are  explained  in  reference  3. 
The  values  of  these  factors  used  in  the  calculation  of  these  fuses  are  kj  =2,  k2  =1, 
and  a  =  5.9  x  10'4  ,  and  were  obtained  from  Di  Marco's  work  (reference  4). 

A  different  reference  size  was  calculated  for  each  fuse  inductance,  since  the 
Maisonnier  size  is  a  function  of  inductance.  Nominal  values  for  bank  capacitance 
and  foil  thickness  were  used  to  calculate  the  reference  sizes.  The  actual  values 
were  subsequently  measured,  and  the  reference  sizes  used  are  about  15%  smaller 
than  they  should  have  been.  The  length  of  the  fuses  varied  from  26.9  to  70.9  cm, 
and  the  width  varied  from  15.0  to  39.4  cm.  The  inductance  of  the  fuses  varied  from 
5. 8  to  37. 9  nH;  however,  most  of  the  fuses  had  an  inductance  of  about  14  nH. 

The  fuses  subdivide  into  two  classes:  those  where  the  ratio  of  the  length  to  the 
width  was  constant,  and  those  where  the  width  was  increased  by  the  same  percentage 
that  the  length  was  decreased.  The  former  group  had  a  nearly  constant  fuse  induct¬ 
ance,  while  the  latter  maintained  a  roughly  constant  total  fuse  material. 

A  coil  of  the  type  normally  used  for  Rogowski  coils  was  used  to  measure  the 
dl/dt.  The  coil  was  placed  in  a  groove  in  the  ground  transmission  plate  and  extended 
several  inches  past  the  edge  of  the  plate.  The  coil  was  not  closed  as  a  Rogowski 
coil  would  be,  but  the  field  away  from  the  parallel  plate  transmission  line  falls  off 
rapidly,  and  the  signal  from  the  coil  should  be  nearly  the  same  as  that  from  a  true 
Rogowski  coil.  In  any  case,  the  coil  sensitivity  was  measured  by  firing  short  circuit 
shots  and  recording  the  coil  output. 

A  passive  RC  integrator  was  used  to  integrate  the  output  of  the  dl/dt  sensing 
coil,  and  the  resulting  signal  was  recorded  on  an  oscillogram.  The  integrator 
employed  a  coaxial  capacitor  to  give  it  good  high-frequency  response.  The  time 
constant  of  the  integrator  is  117  pis,  while  the  foil  always  fused  within  7  pts  after 
the  bank  was  fired. 
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Fuse  voltage  was  measui'ed  by  placing  a  copper  sulfate  water  resistor  across 
the  fuse,  and  measuring  the  current  through  the  resistor  with  a  Tektronix  CT-2  cur¬ 
rent  transformer.  This  was  taken  to  be  proportional  to  the  voltage  across  the  resistor. 
The  resistance  of  the  water  resistor  was  measured  in  two  ways.  A  low  voltage  1-lcHz 
AC  impedance  bridge  (General  Radio  model  1603)  was  used.  Also,  a  200-V,  1-fis 
pulse  from  a2D-21  thyratron  discharge  was  applied  across  the  resistor.  When  the 
resistor  was  new,  these  two  agreed  to  within  the  limits  of  the  measurements  (~5%). 

As  the  resistor  aged,  it  developed  bubbles,  which  were  replaced  with  deionized  water. 
This  did  not  seem  to  effect  the  resistance  very  much  as  measured  using  the  fast  pulses. 
Over  a  period  of  about  1  year,  the  resistance  increased  by  about  10%.  However,  the 
low  frequency  resistance  (1-kHz)  increased  by  over  100%.  Therefore,  the  high 
frequency  resistance  was  measured  regularly  during  the  experiments.  Moreover, 
the  first  fuse  size  that  was  fired  in  a  given  set  of  experiments  was  repeated  at  the 
end  of  the  experiments.  No  significant  difference  was  detected  in  the  voltage  magni¬ 
tude. 

The  frequency  response  of  the  voltage  probe  can  be  assessed  by  observing  the 
rise  of  the  fast  discharge  used  to  calibrate  it.  The  risetime  of  the  pulse  from  the 
thyratron  tube  was  less  than  30  ns,  and  the  risetime  of  the  probe  output  pulse  was 
also  less  than  30  ns.  No  attempt  was  made  to  measure  the  frequency  response  more 
accurately. 

Both  current  and  voltage  were  recorded  as  functions  of  time  on  oscilloscopes. 

In  addition,  one  oscilloscope  had  the  current  signal  applied  to  the  vertical  input,  and 
the  voltage  signal  applied  to  the  external  horizontal  input.  Delays  were  introduced 
into  the  signals  to  ensure  that  simultaneous  signals  on  the  input  cables  reached  their 
respective  deflection  plates  on  the  cathode-ray  tube  simultaneously  (within  2  ns). 

This  gave  a  convenient  way  to  define  the  relative  positions  of  the  current  and  voltage 
pulses.  Since  the  signals  are  changing  very  rapidly  when  the  foil  fuses,  it  is  impor¬ 
tant  to  know  the  relative  timing  accurately. 
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EXPERIMENTAL  RESULTS 

Figure  3  shows  oscillograms  of  the  data  taken  when  a  nominal  Maisonnier-size 
fuse  (46.6-cm-long  by  27.1-cm-wide)  was  fired.  The  top  trace  is  the  voltage,  and 
the  bottom  is  the  current.  The  peak  voltage  of  about  210  kV  occurred  at  about  4.4  us 
into  the  discharge.  A  maximum  of  about  1.7  MA  was  flowing  through  the  fuse  at  3.4 
Us. 

The  data  were  digitized  by  reading  the  oscillograms  on  a  digitizing  machine. 
Primary  analysis  was  then  carried  out  using  the  CDC  6600  computers  at  the  Air 
Force  Weapons  Laboratory. 

Figures  4  and  5  show  voltage  and  current  data  for  representative  fuses  of  the 
constant  inductance  family.  While  the  inductance  changes  only  from  13  to  15  nH,  the 
mass  of  copper  changes  from  9.7  to  53.0  g.  Since  the  inductance  is  nearly  the  same, 
it  is  not  surprising  that  the  initial  voltage  is  nearly  the  same  for  all  the  fuses.  More¬ 
over,  the  current  is  nearly  the  same  until  the  fuse  resistance  increases,  causing  the 
current  to  shut  off.  It  is  also  not  surprising  that  the  larger  mass  foils  fuse  1;  ter  and 
at  higher  currents. 

Figures  6  and  7  show  the  performance  of  fuses  that  have  nearly  the  same  foil 
mass  but  very  different  inductances.  The  Inductances  vary  from  6  to  38  nH,  but  the 
mass  varies  only  from  22.7  to  27.0  g.  The  nominal  Maisonnier  size  fuse  (46.6-cm- 
long  by  27.1-cm-wide)  has  the  largest  mass,  and  the  two  extreme  foils  (26.9-cm- 
long  by  39.4-cm-wide  and  70.9-cm-long  by  15.0-cm-wide)  have  the  smallest  mass. 

It  is  clear  that  the  inductance  of  the  fuse  plays  a  significant  roil  in  the  initial  fuse 
voltage  and  bank  dl/dt.  For  these  fuses,  the  smaller  the  cross-sectional  area,  the 
earlier  it  fuses. 

Both  families  shown  the  peak  voltage  decreasing  as  the  cross-sectional  area 
of  the  fuse  increases,  except,  in  each  case,  for  the  narrowest  fuse.  Figure  8 
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i-Long,  27.1-Cm-Wide  Copper  Foil 
discharge.  Peak  Current  of  1 . 66 


VOLTAGE  (kV) 


Figure  4.  Voltage  Data  for  Approximately  14-nH  Copper  Foil  Fuses. 

(Solid:  27.9  cm  long  by  16.3  cm  wide;  dashed:  37.2  cm 
long  by  21.7  cm  wide;  long  dash:  46.6  cm  long  by  27.1 
cm  wide;  dot -dash:  55.9  cm  long  by  32.5  cm  wide;  and 
double  dot -dash:  65.2  cm  long  by  38.0  cm  wide.) 
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Figure  5.  Current  Data  for  Approximately  14-nH  Copper  Foil  Fuses. 
(Same  sizes  as  for  figure  4.) 
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Figure  6.  Voltage  Data  for  Approximately  25-g  Copper  Foil  Fuses. 

{Solid:  70.9  cm  long  by  15.0  cm  wide;  dashed:  57.7  cm 
long  by  21.0  cm  wide;  long  dash:  4.6  cm  long  by  27.1 
cm  wide;  dot -dash:  36.5  cm  long  by  33.2  cm  wide;  and 
double  dot -dash:  26.9  cm  long  by  39.4  cm  wide.) 


PEAK  VOLTAGE  (kV) 


CROSS  SECTION  (cm2) 


Figure  8.  Peak  Fuse  Voltage  Plotted  Against  Foil  Cross  Section. 
(  a  :  constant  foil  mass  fuses 
a:  constant  inductance  fuses 
The  line  is  a  cubic  least-squares  fit  to  all  the  points.) 
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shows  the  peak  voltage  plotted  against  the  cross-sectional  area.  The  solid  line  is  a 
cubic  least-squares  fit  to  the  points.  Clearly,  a  single  curve  can  reasonably  fit  the 
data  for  both  families  of  fuses.  A  Maisonnier-size  fuse,  with  an  inductance  of  14  nH 
for  the  bank  configuration  used  here,  would  have  a  cross-sectional  area  of  0.078  cm2. 

Data  in  figures  9  and  10  are  the  voltage  and  current  traces  for  two  fuses  of 
equal  length  and  cross-sectional  area,  and,  therefore,  mass  of  copper.  The  differ¬ 
ence  is  that  one  was  split  into  halves  which  were  separated  in  space.  A  46.6-cm- 
long  by  27.1-cm-wide  foil  was  split  lengthwise  into  two  equal  pieces.  These  were 
then  attached  to  the  fuse  clamp  on  the  bank  in  the  usual  way.  The  inside  edges  of  the 
two  pieces  were  separated  by  18.7  cm.  This  set  of  two  parallel  fuses  was  then  fired. 

A  comparison  of  the  data  in  figures  9  and  10  show  that  the  split  foil  fused  sub¬ 
stantially  earlier,  and  reached  significantly  higher  peak  voltage  and  current.  The 
difference  in  inductance  between  the  split  and  normal  fuses  is  calculated  to  be  slightly 
less  than  1  nH,  which  is  less  than  2%  of  the  total  system  inductance.  Therefore,  the 
inductance  difference  alone  cannot  explain  the  differerence  in  performance. 

A  further  test  was  made  by  splitting  a  fuse  and  then  putting  the  two  halves 
together  with  no  space  between  them.  This  resulted  in  a  shot  that  was  not  signifi¬ 
cantly  different  from  the  single-piece  fuses  of  the  same  size.  These  tests  are  taken 
to  indicate  that  edge  effects  are  very  important  to  the  fuse  behavior,  which  is  not 
surprising  since  the  current  density  is  known  to  be  higher  at  the  edges  of  parallel 
plate  transmission  lines  than  at  the  center. 

Voltage  oscillations  on  the  early  time  data  for  the  split  fuse  were  unique  to 
that  particular  configuration.  However,  they  were  present  on  all  of  the  two-piece 
fuses  that  were  fired.  These  oscillations  were  not  present  on  the  three-  and  four- 
piece  fuses  that  were  fired. 

Figures  11  and  12  compare  the  two-piece  and  four-piece  fuses.  Total  cross- 
sectional  area  and  length  of  the  foils  were  the  same:  0.065  cm2  and  46.6  cm, 
respectively.  Aside  from  the  early  time  oscillations,  there  is  no  significant  differ¬ 
ence  in  fuse  behavior. 
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Figure  9. 


Comparison  of  the  Voltage  Data  for  One-Piece  and  Two-Piece 
Fuses.  (Solid  line:  one-piece  fuse;  dashed  line:  two-piece 
fuse.  Both  fuses  were  46.6  cm  long  and  have  a  cross-sectional 
area  of  0.065  cm5.) 


.00 


.70  1  .40  2.  10  2. BO  3.50  4.20  4.90  5.60  6.30  7.00 

TIME  (/i sbc) 

Figure  10.  Comparison  of  the  Current  Data  for  One-Piece  and  Two- 
Piece  Fuses.  (Rest  of  caption  same  as  for  figure  9.) 
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Figure  11.  Comparison  of  the  Voltage  Data  for  Two-Piece  and  Four- 
Piece  Fuses.  (Solid  line:  two-piece  fuse;  dashed  line:  four 
piece  fuse.  Both  fuses  were  46.6  cm  long  and  had  a  cross 
sectional  area  of  0.065  cm*.) 


There  was  6. 3-cm  separation  between  each  of  the  fuse  elements  in  the  four- 
piece  fuse.  Consequently,  both  fuse  packages  had  an  overall  width  of  about  46  cm. 
Perhaps  this  is  the  important  parameter;  nevertheless,  it  was  not  investigated. 

In  order  to  check  the  overall  reproducibility  of  the  data,  the  last  fuse  fired  was 
the  same  size  as  the  first  one.  Figures  13  and  14  compare  the  results  of  these  two 
fuses.  There  is  a  significant  difference;  nevertheless,  the  differences  are  much  less 
than  those  between  different  fuse  sizes  (c.p.  figures  4-7).  Even  splitting  the  fuse  is 
much  more  significant.  It  is  felt  that  the  differences  arise  from  differences  between 
the  copper  at  the  beginning  and  the  end  of  the  roll. 

Instantaneous  fuse  resistance  can  be  calculated  using  the  equation 


HF(t> 


VF(t)  -  LrI(t) 
I(t) 


(1) 


where  V  is  the  fuse  voltage,  L_  is  the  fuse  inductance,  I  is  the  current,  and  I  is 
x  r 

the  time  rate  of  change  of  the  current.  A  quantity  that  has  units  of  resistivity  can 
be  defined  by 


s  RF(t>  f 


(2) 


where  s  is  the  initial  cross-sectional  area  of  the  fuse  and  £  is  the  length.  Clearly, 
P  is  the  average  resistivity  of  the  fuse  material  if  the  foil  dimensions  have  not 
changed.  Henceforth,  P  will  be  treated  as  though  It  were  the  resistivity  of  the 
copper. 
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Figure  13.  Comparison  of  the  Voltage  Data  for  Fuses  Fired  Early 
(Solid  Line)  and  Late  (Dashed  Line)  in  the  Experimental 
Program.  (Both  fuses  were  46.6  cm  long  and  27.1  cm 
wide.  Approximately  5  months  and  50  shots  separated 
these  two. ) 
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Resistivity  is  plotted  as  a  function  of  specific  energy  absorbed  by  the  fuse  mate¬ 
rial  in  figure  15  for  five  different  fuses  of  the  constant  inductance  family.  Specific 
energy  is  defined  by 


where 


<?(t) 


t 


P(t')  J 2(t')  dt' 


(3) 


J(t) 


(4) 


and  y  is  the  mass  density  of  copper  (8. 96  g/cm3 ,  reference  9). 

Several  features  of  these  curves  are  worth  noting.  As  the  peak  current  density 
increases,  the  specific  energy  at  which  the  resistivity  starts  to  climb  steeply  increases. 
Moreover,  the  climb  is  steeper  for  lower  peak  current  density.  Particularly  for  high 
current -density  fuses,  there  is  evidence  of  a  plateau  in  the  resistivity  versus  specific 
energy  curve.  This  plateau  is  at  a  higher  resisitivity  for  lower  current  densities. 
Similar  results  have  been  reported  for  exploding  copper  wires  (reference  10). 

Resistivity  of  the  same  fuses  is  plotted  against  specific  action  in  figure  16. 

Specific  action  is  defined  by 

t 

G(t)  =  j  Jit')1  dt'  (5) 

o 


The  specific  action  at  which  the  resistivity  climbs  steeply  is  seen  to  be  a  function  of 
current  density.  Nevertheless,  the  dependence  is  somewhat  weaker. 


RESISTIVITY  (*/ohm-cm) 


Figure  15.  Resistivity  Plotted  Against  Specific  Energy  Absorbed  for  Five 
Fuses  of  the  Constant  Inductance  Family.  (Solid:  27.9  cm  long 
by  16.3  cm  wide;  dashed:  37.2  cm  long  by  21.7  cm  wide; 
long  dash:  46.6  cm  long  by  27.1  cm  wide;  dot -dash:  55.9  cm 
long  by  32.5  cm  wide;  and  double  dot-dash:  65.2  cm  long  by 
38.0  cm  wide.) 
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RESISTIVITY  (^ohm-cm) 


Figure  16.  Resistivity  Plotted  Against  Specific  Action  for  Five  Fuses  of  the 
Constant  Inductance  Family .  (Rest  of  caption  same  as  for 
figure  15.) 
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SECTION  IV 


RESISTIVITY  MODEL 

In  general  the  resistivity  versus  energy  density  curves  for  the  fuses  within  each 
family  can  be  ordered  by  the  peak  current  density.  That  is,  the  foil  that  reaches  a 
higher  peak  current  density  fuses  at  a  higher  specific  energy.  If  fuses  from  the  con¬ 
stant  mass  family  are  compared  to  the  fuses  from  the  constant  inductance  family, 
this  ordering  is  no  longer  possible. 

Figure  17  shows  the  resistivity  versus  specific  energy  curves  for  three  fuses. 
Of  the  three  fuses,  one  belonged  to  the  constant  inductance  family,  one  belonged  to 
the  constant  mass  family,  and  one  belonged  to  both  families.  These  three  fuses 
reached  peak  densities  of  30.3,  28.3,  and  25.6  MA/cmJ  ,  respectively.  Clearly, 
these  curves  cannot  be  ordered  on  the  basis  of  peak  current  density. 

Figure  18  shows  the  resistivity  versus  specific  action  curves  for  the  same 
three  fuses.  These  curves  can  be  ordered  on  the  basis  of  peak  current  density. 

When  all  of  the  data  are  compared,  the  ordering  is  no  longer  perfect.  Nevertheless, 
it  is  much  better  than  trying  to  order  the  resistivity  versus  specific  energy  curves. 

It  is  tempting  to  try  to  find  a  simple  dependence  of  the  burst  specific  action  on 

current  density.  (Burst  specific  action  is  defined  as  that  specific  action  at  which  the 

copper  has  completely  vaporized.)  In  other  words,  is  there  a  parameter  of  the  fused 

circuit  behavior,  G  ,  such  that  resistivity  as  a  function  of  G  can  be  uniquely  defined 
c  C 

without  an  explicit  dependence  on  J  ?  We  call  Gg  the  effective  specific  action  and 
define  it  like  specific  action  with  the  addition  of  an  implicit  functional  dependence  on 
J  called  F(J) 

t 

°e(t>  Uf  j2{V)  F(J(t,))  dt'  (6) 

o 
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RESISTIVITY  (//ohm -cm) 


Figure  17.  Resistivity  Plotted  Against  Specific  Energy  for  Three  Fuses. 

(Solid  line:  37.2  cm  long  by  21.7  cm  wide,  14  nH,  17.3  g  of 
copper,  30.3  MA/cmJ  peak  current  density,  dashed  line:  57.7 
cm  long  by  21.0  cm  wide,  23  nH,  25.9  g  of  copper,  23. 3  MA/cm2 
peak  current  density;  dot -dash:  46.6  cm  long  by  27.1  cm  wide, 

15  nH,  27  g  of  copper,  25.6  MA/cm2  peak  current  density.) 
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The  function  F(J)  is  a  correction  factor  that  must  be  employed  when  using  some 
resistivity  fur  ion.  T.  J.  Tucker  and  R.  P.  Toth  have  published  resistivity  versus 
specific  action  curves  for  many  materials  (reference  11).  These  data  were  taken  at 
constant  current  densities,  and  they  are  ideal  as  a  basis  for  a  resistivity  model. 

The  problem  now  is  to  find  the  function  F(J),  such  as  that  at  the  time  of  peak 
voltage, 


G  /t\  =  C  =  constant  (7) 

e\  B  /  o 


for  all  the  fuses.  Several  different  functional  forms  were  tried.  The  form  that  was 
both  simple  and  did  a  good  job  in  fitting  the  data  is  the  linear  form 


F(J)  ■  1  +  A^J  - 


(8) 


A  and  J  are  to  be  taken  from  fits  to  the  data, 
o 

Tucker  reports  the  peak  in  the  resistivity  curve  for  copper  to  occur  at  1730 
MA2  -//s/cm4 .  Using  this  to  be  the  effective  specific  action  at  which  copper  bursts, 
a  fit  to  the  data  presented  here  leads  to  the  values  for  A  and  Jq  shown  in  table  1. 

TABLE  1 

RESISTIVITY  MODEL  CONSTANTS 

1730  MA2  -//s/cm4 
-0.0236  cm2 /MA 
25.7  MA/cm* 


G  * 
o 

A 

J 

o 


♦Reference  11. 
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The  resistivity  versus  specific  action  curve  is  steeper  when  the  foils  fuse  at 
lower  specific  actions.  Exactly  how  to  do  this  with  a  model  is  a  problem,  if  the  model 
consists  of  a  data  table .  Therefore,  it  was  decided  to  fit  a  curve  to  Tucker's  data. 
The  curve  that  was  chosen  consists  of  the  sum  of  eight  gaussian  functions.  The  cen¬ 
troids,  heights,  and  full  widths  at  half  maximum  (FWHM)  are  listed  in  table  2. 


TABLE  2 

PARAMETERS  FOR  THE  FIT  TO  TUCKER'S  DATA* 

Center  Height  FWHM 

(MA2  -j/s/ cm4 )  (uQ-cm)  (MA2  -us /cm*  ) 


6589 

3707 

3994 

1733 

516.1 

66.9 

966 

4.8 

93.6 

1084 

2.9 

181.2 

1588 

8.5 

73.2 

1594 

10.2 

28.6 

1825 

56.6 

116.9 

♦Reference  11. 

As  the  foil  fuses  at  smaller  values  of  specific  action,  the  rise  becomes  steeper, 
and  the  plateau  in  the  resistivity  occurs  at  higher  values.  This  can  be  accomplished 
with  a  gaussian  by  increasing  the  amplitude  and  decreasing  the  FWHM.  It  was  decided 
to  change  the  amplitude  and  FWHM  by  the  same  ratio;  that  is,  the  area  under  each 
gaussian  in  the  fit  was  kept  constant.  A  reasonable  fit  to  the  data  was  obtained  by 
miltiplying  the  amplitude  of  each  gaussian  by  the  factor 


1 


+ 


1.2 


(‘ 


(9) 
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where  G  is  the  specific  action  and  Gg  is  the  effective  specific  action.  The  FWHM  are 
then  divided  by  this  factor. 

Copper  vapor  from  the  fuse  quickly  condenses  on  the  glass  beads  inside  the  fuse 
package.  Therefore,  the  fuse  normally  does  not  restrike,  and  the  resistance  remains 
high.  In  order  to  model  this,  the  resistivity  is  kept  constant  after  the  effective  specific 
action  has  reached  1734  MAJ-/is/cm‘. 

Predicted  resistivity  as  a  function  of  specific  action  is  plotted  for  three  current 
densities  in  figure  19.  These  curves  were  made  using  the  model  described  above, 
assuming  cons'. ant  current  densities.  Since  the  highest  current  density  obtained  in 
the  present  data  was  about  35  MA/cm2,  it  is  unwise  to  extrapolate  this  resistivity 
model  much  beyond  that  level.  The  resistivity  curves  in  figure  19  show  regions  where 
the  resistivity  is  decreasing.  This  is  an  artifact  of  the  way  in  which  the  shape  of  the 
resistivity  curve  is  modified  and  it  is  a  shortcoming  of  this  resistivity  model. 

In  order  to  see  how  well  the  use  of  the  effective  specific  action  as  the  independ¬ 
ent  variable  removes  the  current  density  dependence  from  the  resistivity,  compare 
figures  20  and  21.  Resistivity  is  plotted  as  a  function  of  specific  action  in  figure  20 
for  three  fuse  sizes.  Resistivity  is  plotted  as  a  function  of  effective  specific  action 
for  the  same  three  fuses  in  figure  21 .  The  upper  end  of  each  curve  corresponds  to 
peak  voltage  in  that  fuse. 
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Figure  19.  Model  Predictions  for  Resistivity  Versus  Specific  Action  at 

Three  Current  Densities.  (Solid:  30  MA/cm1;  dash:  20  MA/cm* 
dot-dash:  10  MA/cm1.) 
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Figure  20.  Plot  of  Resistivity  Versus  Specific  Action  for  Three 
Fuses.  (Solid  line:  27.9  cm  long  by  16.3  cm  wide; 
dashed  line:  31.7  cm  long  by  36.3  cm  wide;  dot- 
dashed  line:  41.9  cm  long  by  24.4  cm  wide.) 
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SECTION  V 


FUSE  CODE  AND  COMPARISONS  WITH  EXPERIMENT 

A  computer  code  was  written  to  simulate  the  performance  of  copper  foil  fuses. 
The  resistivity  model  described  in  the  previous  section  is  built  into  the  code.  This 
code  has  been  used  to  predict  the  behavior  of  a  wide  range  of  fuse  sizes  with  generally 
good  results. 

Since  edge  effects  are  important  to  fuse  behavior,  some  two-dimensional  effects 
are  explicitly  treated.  In  the  code,  the  fuse  is  considered  to  be  a  number  of  parallel 
copper  ribbons  with  the  edge  of  one  touching  the  edge  of  the  next,  and  all  connected  in 
parallel  to  the  bank.  Typically,  the  fuse  is  thought  of  as  consisting  of  100  ribbons. 
Mutual  inductances  are  calculated  and  the  ribbons  interact  magnetically,  but  current 
is  not  allowed  to  transfer  from  within  one  ribbon  to  another. 

It  is  assumed  that  there  is  the  same  voltage  drop  around  each  loop  containing  a 
ribbon.  At  any  given  time,  therefore,  this  voltage  is  given  by 


e  =  Vi  +  Vi  +J3  “iA  <10> 

where  R.,  L.,  I.,  and  i  are  the  resistance,  self  inductance,  current,  and  time  rate 
of  charge  of  current  in  the  ith  ribbon,  is  the  mutual  inductance  between  ribbons 
i  and  k,  and  is  the  time  rate  of  change  of  current  in  the  kth  ribbon.  If  we  define 
=  M..,  then  this  equation  becomes 


*m*lh 


(U) 
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There  are,  of  course,  N  equations,  where  N  is  the  total  number  of  loops.  The 
external  circuit  equation  is 


~  Q  =  R  I  +  L  I  +  £ 
C  x  x 


(12) 


where  C  is  the  bank  capacitance,  Q  is  the  charge  on  the  bank,  Rx  and  are  the 

external  circuit  resistance  and  inductance,  respectively,  I  is  the  total  current  and 
• 

I  is  the  time  rate  of  change  of  total  current.  Combining  these  two  equations  results 
in  a  set  of  N  simultaneous  linear  differential  equations  that  must  be  solved. 

These  simultaneous  equations  compose  an  Nth  order  matrix  equation,  which  is 

solved  for  the  I  at  a  succession  of  time  steps.  At  each  time  step  the  Q,  I,  and  I. 

^  • 

are  calculated  from  the  previous  values  of  Q,  I,  I.,  and  I..  Then  the  present  values 
of  I.  are  calculated.  The  resistivity  model  is  used  to  calculate  the  R.  at  each  time 
step.  This  procedure  is  repeated  until  the  problem  is  finished.  Time  step  size  is 
varied  to  ensure  that  I  changes  only  by  a  small  amount  in  each  time  step.  More 
detailed  information  on  the  code  can  be  found  in  reference  12. 

Simulations  were  run  for  all  of  the  foil  sizes  used,  including  split  foils. 

Figure  22  is  a  comparison  between  the  simulations  and  experiments  for  four  different 
fuse  sizes.  It  can  be  seen  that  the  simulations  range  from  rather  good  to  only  fair. 

In  all  cases,  the  initial  voltage  of  the  simulation  is  lower  than  the  experiment. 
However,  this  difference  does  not  have  the  same  ratio  for  all  the  foil  sizes.  There¬ 
fore,  it  is  not  the  result  of  a  simple  calibration  error  in  the  voltage  probe. 

A  curve  was  least  square  fit  to  the  dl/dt  data  of  each  fuse  shot.  Inductance  was 
calculated  from  the  intercept  of  each  curve  and  the  calculated  fuse  inductance  sub¬ 
tracted.  The  bank  inductance  was  calculated  by  averaging  these  numbers  to  be  44 
+  3  nH.  Therefore,  fuse  inductance  calculations  are  not  in  error  by  enough  to 
account  for  this  discrepancy. 
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Figure  23  compares  fuse  voltage  curves  near  peak  voltage  for  two  simulations. 
The  solid  curve  is  for  a  normal  one-piece  fuse,  while  the  dashed  line  is  for  a  split 
fuse.  There  is  a  difference  between  the  two  curves;  nevertheless,  it  is  very  small 
compared  to  the  experimental  difference  (c.p.  figure  9). 

Although  the  fuse  code  does  predict  qualitative  fuse  behavior,  it  does  not  make 
predictions  that  are  as  accurate  as  one  would  like.  It  is  felt  that  the  reason  for  this 
is  that  two-dimensional  effects  are  very  important,  that  is,  significant  current  can 
flow  perpendicular  to  the  length  of  the  fuse.  The  present  fuse  code  does  not  allow 
for  these  currents.  Perhaps  a  code  such  as  PLATE  (reference  13)  or  EBFI  (reference 
14)  with  the  present  resistivity  model  would  be  more  successful. 


VOLTAGE  (kV) 


Figure  23.  Comparison  of  Computer  Simulation  Voltage  Pulses  for  One-Piece 
and  Two-Piece  Copper  Foil  Fuses.  (Both  foils  are  4G.6  cm  long 
with  a  cross-sectional  area  of  0.065  cm2.) 
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SECTION  VI 


CHOICE  OF  FUSE  SIZE 

In  order  to  select  the  fuse  size  that  is  appropriate  for  a  given  capacitor  bank, 
we  adopt  the  first  of  Maisonnier's  criteria.  That  is,  the  fuse  should  burst  at  peak 
current  (reference  3).  Furthermore,  it  is  assumed  that  before  the  foil  bursts,  its 
resistance  is  negligible  and  that  the  current  density  is  uniform  everywhere  in  foil. 

Under  these  circumstances,  the  current  waveform  will  be 


I  =  V  C  cj  sin  U it 
o 


(13) 


where  C  and  are  the  bank  capacitance  and  initial  charge  voltage,  and 


U  = (LC) 


-1/2 


(14) 


where  L  is  the  total  inductance.  The  condition  for  burst  at  peak  current  is 


G 

o 


G 

e 


n/  20) 

•/  j,(i+a 


o 


(J 


(15) 
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where  Gq,  A,  and  Jq  are  const;ints  and  J  is  the  current  density.  Substituting  i/s  for 
J  and  integrating,  gives 


o  (  o)  U L  4  s3  L  3  s3 


Multiplying  by  7; —  and  collecting  the  terms  give 


1-AJ  V3C  „  V3C2 

s3 _ o_  7T  o  2  _A_  __o _ 

G  4  0/L  3  G  L 

o  0 


If  there  were  no  dependence  of  the  burst  specific  action  on  the  current  density,  then 
A  would  be  0.  Solving  equation  (17)  for  s  with  A  =  0  gives 


4Go  (LC|lA 


Maisonnier's  equation  for  the  cross-sectional  area  of  the  fuse  can  be  written  in  the 
form 


2  k,  a 


where  k  and  a  are  parameters  that  are  defined  in  reference  3.  By  comparing 
equations  (18)  and  (19),  we  can  make  the  identification 


i 


a 


4G 


(20) 


Using  the  value  5.9  x  1014  for  a  (reference  4),  gives  kj  =  1.4.  Di  Marco  and 
Burkhardt  found  that  a  value  of  about  2  gave  a  good  fit  to  their  data  (reference  4). 

The  present  data  indicates  that  the  burst  specific  action  is  a  function  of 
current  density  and  A  is  not  0.  Defining 


and 


1  -  AJ 


a  =  _ 


7 T 

4 


V2C 

o 

wL 


b 


2  A_ 

3  G 

o 


(21) 


(22) 


the  cubic  equation 


s3  +  as  +  b  =  0 


(23) 


results. 


If  b2/4  +  a  3/27  >  0,  there  will  be  only  one  real  solution  to  equation  (23). 
Conversely,  if  it  is  less  than  0,  there  will  be  three  real  roots.  The  dividing  line 
between  these  two  cases  occurs  when 


0 


A2V6C4 

o 

9G2L2 

o 


(l  -  AjJV3VoC 

172SG3<oL2 

o 


(24) 


This  can  be  reduced  to 

ttVi -AJ  V 

W  =  192G  A2  '  =  6*95xl°5  s"  (25> 

o 

This  corresponds  to  a  quarter-cycle  time  of  2.26  fis. 

If  the  natural  risetime  of  the  bank  is  slower  than  this,  there  will  be  three  real 
solutions.  They  are 


where  n  =  0, 1, 2.  The  angle  <p  is  defined  by 


' 1 


192G  A2 


)■ 


cos  <t>  =  - 


(27) 


.1 


Clearly,  Q  will  always  be  greater  than  tt/2  and  less  than  n.  Therefore,  ~  will 
always  be  greater  than  7r/6  and  less  than  7r/3,  and  there  will  always  be  two  positive 
real  solutions  for  s.  For  example,  if  L  =  60  nH  and  C  -  158  MF ,  0  =  2.32  rad,  and 
the  three  solutions  for  s  are  0.100  cm2 ,  -0.135  cm1,  and  0.035  cm2.  A  comparison 
with  figure  8  shows  that  the  two  positive  roots  bracket  the  experimental  fuses. 

If  the  assumption  of  zero  resistance  until  the  foil  fuses  were  strictly  valid,  the 

fuses  with  cross-sectional  areas  0.100  and  0.035  cm2  would  reach  peak  current 

densities  of  25.7  and  73.5  MA/cm2,  respectively.  As  presented  here,  F(J)  becomes 

negative  above  68  MA/cm2.  This  means  that  for  the  smaller  positive  root,  G  (t)  is 

not  monotonic.  In  fact,  G  ft)  mcr eases  to  values  larger  than  G  and  then  decreases 

e  o 

until  its  final  value  is  G  .  This  is  clearly  unphysical,  since  the  foil  would  fuse  the 

o 

first  time  G  (t)  reached  G  .  Indeed  this  is  what  happens  experimentally.  Consequently 
e  o 

only  one  of  the  roots  is  physically  meaningful. 

If  the  risetime  of  the  bank  is  faster  than  2. 26  //s,  there  is  only  one  real  root  to 
equation  (23) ,  and  it  is 


Since  A/Gq  is  negative,  s  will  always  be  negative. 


(28) 
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Clearly,  there  is  little  meaning  to  the  solution  of  equation  (23)  for  fast  rise 
banks. 

In  order  for  the  correction  factor  F(J)  to  be  valid,  J  must  not  greatly  exceed 
the  maximum  value  of  J  that  was  encountered  in  the  experiments  from  which  F(J) 
was  inferred.  That  is,  J  must  not  greatly  exceed  35  MA/cm2.  Equations  (13)  and 
(26)  can  be  used  to  infer  the  peak  current  density 


I 


3G 


3EZ 


AJ 


4  1  <t> 

LC  cos  f 

u 


(29) 


If  the  risetime  is  3  fxs. 


J  .  -  36  MA/cm2 
pk 


(30) 


Moreover,  if  the  risetime  is  longer  than  3  fjs,  the  peak  current  density  will  be 
smaller  than  36  MA/cm2.  Therefore,  the  optimum  fuse  cross-sectional  area  can  be 
predicted  if  the  risetime  is  longer  than  about  3  ^s.  This  can  be  expected  to  be 
somewhat  better  than  Malsonnier's  equation,  because  it  takes  the  current  density 
dependence  of  the  fuse  behavior  into  account. 

In  order  to  select  the  length  of  the  fuse,  one  more  criterion  must  be  specified. 
Since  the  fuse  resistivity  plateaus,  It  is  convenient  to  choose  the  shut-off  time  r. 
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Then  the  current  in  the  fuse  will  decay  in  a  time  characteristic  of  L/R,  where  L  is 
the  total  inductance  and  R  is  the  plateau  resistance  of  the  fuse.  Therefore,  the  length 
is  specified  by 


(31) 


where  pQ  is  the  plateau  resistivity. 

A  reasonable  value  for  the  plateau  resistivity  is  300  /jQ-cm  (figure  15).  If  one 
chooses  r  to  be  300  ns,  L  is  60  nH,  and  s  is  0.065  cm2,  then 


l  =  43  cm 


(32) 


This  is  very  close  to  the  0.065  cm2  by  46.6  cm  fuse  that  was  chosen  using  Maisonnier’s 
energy  dissipation  criterion  (reference  3).  Notice,  however,  that  one  can  choose  a 
wide  range  of  lengths  to  select  the  desired  shut-off  time  without  regard  to  the  amount 
of  energy  being  absorbed. 
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SECTION  VII 


CONCLUSION 


Copper  foils  have  been  successfully  used  as  fuses  on  a  200-kj  capacitor  bank. 
While  no  parallel  load  was  switched  into  the  circuit,  the  current  in  the  fuse  was  shut 
off  at  rates  as  high  as  4  x  10,J  A/s.  Voltage  pulses  with  a  full  width  at  half  maximum 
of  about  300  ns  and  peaks  at  4  /is  into  the  discharge  have  been  observed.  Therefore, 
ratios  of  conduction  time  to  turn-off  time  in  excess  of  10  have  been  observed. 

Resistivity  as  a  function  of  specific  action  and  current  density  has  been  inferred. 
In  general,  a  higher  peak  current  density  results  in  fusing  at  a  higher  specific  action. 
The  current  density  dependence  has  the  effect  of  extending  the  range  of  foil  sizes  which 
will  work  as  fuses. 

Postulating  that  the  foil  should  reach  burst  specific  action  at  peak  current  leads 
to  a  method  of  determining  the  optimum  cross-sectional  area.  If  the  burst  specific 
action  is  assumed  to  be  independent  of  the  current  density,  the  resulting  relation  is 
the  same  as  that  derived  by  Maisonnier.  On  the  other  hand,  if  the  current  density 
dependence  inferred  from  the  present  data  is  employed,  a  different  equation  is 
applicable  fuse  size.  This  inferred  optimum  size  is  valid  so  long  as  the  peak  current 
density  does  not  greatly  exceed  35  MA/cm2 .  This  is  equivalent  to  specifying  that  the 
quarter  cycle  risetime  of  the  bank  is  3  /is  or  greater .  Since  the  current  density 
dependence  was  included  in  the  derivation  of  the  present  technique  for  selecting  the 
fuse  cross-sectional  area,  it  should  be  more  appropriate  in  the  region  where  it  is 
valid  than  that  of  Maisonnier. 

Choosing  the  optimum  fuse  length  is  largely  subjective.  The  data  presented 
here  indicate  that  the  fuse  material  will  absorb  large  amounts  of  energy  without 
restriking.  Moreover,  the  foil  will  still  effectively  fuse  at  relatively  low  specific 
energy.  Therefore,  it  is  suggested  that  the  fuse  length  be  chosen  to  give  the  desired 
current  turn-off  time. 
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In  light  of  the  large  range  of  specific  energy  which  a  fuse  will  absorb  and  still 
function,  perhaps  it  is  not  necessary  to  totally  vaporize  the  fuse  material.  Possibly 
all  that  is  necessary  is  to  liquify  the  material  and  drive  it  rapidly  into  the  glass  beads. 
Further  work  would  be  necessary  in  order  to  address  this  question. 
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